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The structure of the flux-tube profile in Abelian-projected (AP) SU(2) gauge theory in the maxi- 
mally Abelian gauge is studied. The connection between the AP flux tube and the classical flux-tube 
solution of the U(l) dual Abelian Higgs (DAH) model is clarified in terms of the path-integral duality 
transformation. This connection suggests that the electric photon and the magnetic monopole parts 
of the Abelian Wilson loop can act as separate sources creating the Coulombic and the solenoidal 
electric field inside a flux tube. The conjecture is confirmed by a lattice simulation which shows 
that the AP flux tube is composed of these two contributions. 
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(NJ ' I. INTRODUCTION 
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CO 

When the QCD vacuum is viewed as a dual superconductor Q, the quark confinement mechanism can be 
3 . immediately understood: the (color-) electric flux associated with a quark-antiquark (q-q) system is squeezed into 
an almost-one-dimensional flux tube by the dual Meissner effect caused by magnetic monopole condensation. This 
picture leads to a linear confinement potential and is a dual analogue of the magnetic Abrikosov vortex in an ordinary 
superconductor 0, 0, S ■ It is natural to expect that it can be quantitatively formulated by a dual version of an 
Abelian Higgs model, the dual Abelian Higgs (DAH) model. The Lagrangian — besides the kinetic terms of each 
field and a minimal coupling between the two fields — should contain a monopole self-interaction term that allows 
for a brokenphase of dual gauge symmetry. The DAH model indeed has an electric flux-tube solution of the static 
q-q system |5j. 

A linear potential emerging from a flux tube is quite welcome to give an interpretation for the area law behavior 
of the Wilson loop observed in lattice QCD simulations Q. It would explain the Regge trajectory pattern or other 
string- like properties of hadrons 0|. Then the problem arises how to derive the dual superconductor scenario from 
QCD, that is, how to formally derive the DAH model from QCD. One also would like to observe certain characteristic 
features of the dual superconductor, such as the formation of flux tubes, through a Monte-Carlo simulation of lattice 
QCD directly. 
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As for the formal derivation, it is known that if magnetic monopoles are introduced as the consequence of Abelian 
projection a la 't Hooft 8] and if the diagonal components of gluons play a dominant role (compared to the off-diagonal 
ones) in the long distance behavior of QCD (Abelian dominance) , a condensed phase of monopoles is realized beyond 
a certain critical scale [H 0|. Remarkably, lattice QCD simulations with non- Abelian configurations undergoing 
't Hooft's Abelian-projection (typically in the maximally Abelian gauge, MAG) support this scenario numerically. 
For instance, the string tension measured by the "Abelian Wilson loop" constructed from the Abelian link variables 
(the "Abelian string tension"), is almost saturating the non- Abelian string tension In this context, applying the 
Zwanziger formalism one can introduce the dual gauge field which is minimally coupled to monopoles. Summing 
over monopole current traj ectories .13, 14] , one can also introduce a monopole field. This formulation finally leads 
to the DAH model \lA lift Il7l Il8j. However, it is difficult to determine the effective couplings of the DAH model 
through this analytical derivation, because one cannot treat the monopole current system quantitatively. In order to 
accomplish this, one would need numerical investigations of monopole dynamics on the lattice, for instance by means 
of the inverse Monte-Carlo method 0, EH 0^ E4 ■ This might require more complicated ansatze for matching the 
monopole actions [2^ . 

Just in order to seek flux-tube configurations in the non-Abelian gauge theory, the profiles of the electric field 
and the monopole current distribution induced by an Abelian Wilson loop have been studied within the Abelian- 
projection scheme 0,EHEHE3- It nas been found that the shapes are similar to those of the flux-tube solution 
in the DAH model. From now we call the former one "Abelian-projected (AP) flux tube" and the latter one "DAH 
flux tube" . We remark that the connection between the AP flux tube and the DAH flux tube is not on equal footing 
because the former contains the quantum effects at work in non-Abelian lattice gauge simulations on the original 
lattice, while the latter is a classical solution obtained by solving the field equations with dual variables. Having in 
mind this conceptual difference, it is still worth to determine the effective couplings of the DAH model, which could 
not be fixed through a formal derivation, through the comparison between the two flux tubes. This is interesting 
because, once the DAH parameters are fixed, one can use the DAH model for further analyses: for discussing hadronic 
objects EH lia. Efl . for investigating the dynamics of the flux tube by deriving an effective string action from the 
DAH model|31iHEHElEHH, etc. 

Up to now, the quantitative status of the comparison between the AP and the DAH flux tubes has not been 
conclusive, although this has been attempted several times 0> EH EH EH EB ■ I n order to find DAH parameters 
which possess physical meaning in this context, at first it is important to understand to what extent the AP flux tube 
can be really related to the DAH flux tube, first of all since they are defined in terms of different (original and dual) 
variables. This should become clear once the duality transformation is carried out in detail. Second, also a more 
systematic study of the AP flux-tube profile is required to have well-controlled lattice data; one needs to check the 
Gribov copy effect hidden in the process of MAG fixing, has to examine to what extent the scaling property is fulfilled, 
should investigate the q-q distance dependence of the flux tube shape etc. on a sufficiently large lattice volume. 

In this paper, we aim to address only the first part, the qualitative and detailed relation between the AP and the 
DAH flux tubes. Here we do not attempt to fix the DAH model parameters. What we plan to do here is to show 
that the AP flux tube has the composed internal structure as the DAH flux tube has, going through the path-integral 
duality transformation of the AP gauge theory. In fact, in the DAH model, as we explain later in detail, the appearance 
of the electric flux tube is due to the superposition of two well distinguished components, a Coulombic electric field, 
directly induced by the electric charges, and a solcnoidal electric field induced by a monopole supercurrent. They 
are responsible for the Coulombic and the linearly rising part, respectively, of the inter-quark potential in the DAH 
model. If the electric flux profile can be uniquely decomposed in the case of the AP flux tube as well, analogously 
to the DAH flux tube, this will be an additional argument in favor of the DAH model description, which will be 
important for further quantitative discussions. 
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The guiding idea to discover this kind of structure also in the AP flux tube comes from the measurement of the q-q 
potential in terms of the Abelian Wilson loop. The investigation of the Abelian Wilson loop using the decomposition 
into an electric photon part ("photon Wilson loop") and a magnetic monop ole p art ("monopole Wilson loop") shows 
that also the Abelian potential consists of a Coulombic and a linear potential |38|, |39|, |4fJ . We notice that this structure 
is quite similar to that of the q-q potential in the DAH model. 

The paper is organized as follows. In section [H] we shall discuss the theoretical connection of the photon and the 
monopole Wilson loops with the composed structure of the DAH flux tube. We do this by closely looking at the 
path-integral duality transformation of the AP gauge theory. Motivated by lattice results on the effective monopole 
action we adopt, as our starting point, a Villain type compact QED as the approximate action of the effective, AP 
gauge theory. In section HTT1 we present the numerical results, the flux profile induced by the photon and monopole 
Wilson loops, measured within SU(2) lattice gauge theory in the MAG. We come to the conclusion that the AP flux 
tube is composed out of Coulomb and solenoidal parts, which add up to the full electric flux tube, in the same manner 
as the DAH flux tube. Section llVl is the summary. 

The due improvement in the systematic study of the AP flux tube including all details of the quantitative analysis 
of our lattice data, along the guiding lines formulated in the present paper, is the subject of our follow-up paper 411. 



II. THE COMPOSED STRUCTURE OF THE FLUX TUBE IN THE DAH MODEL 



In this section, based on a path-integral analysis, we discuss a possible theoretical relation between the electric- 
photon and magnetic-monopole parts of the Abelian Wilson loop in the AP-SU(2) lattice gauge theory and the 
composed internal structure of the flux-tube solution in the U(l) DAH model. 



From lattice studies of the effective monopole action in the MAG [13, |2JJ, |21j, |22j , it is numerically suggested that, 
at some infrared scale, the partition function of the AP-SU(2) theory is represented by the Villain type modification 
of compact QED. Thus, we regard it as the effective AP gauge theory and start from the partition function 
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(F, AD F)+i(9,j) 



(2.1) 



F(C 2 ) is the field strength 



F = d9- 27m (m) 



(2.2) 



which is composed of compact link variables, 9(C\) S [— tt,tt), and magnetic Dirac strings, n^ m \C2) € % 45 
corresponds to the Abelian gauge field, which interacts with an external electric current j(C\) £ The operator D 
is a general differential operator and A is the Laplacian on the lattice. In the infrared limit, it is numerically shown 
that the operator D is well-described by the following form: D = (3 e A~ l +a + ^A, where f3 e , a and 7 are renormalized 
coupling constants of the monopole action which satisfy the relation (3 e ~> a, 7 0|- The (inverse) effective gauge 
coupling is (3 e = 4/e 2 . Since the magnetic Dirac strings n 1 --" 1 ^ are bordered by magnetic monopole currents k{C^) as 
d n ( m ) — (hence dk = 0), the Abelian Bianchi identity is now violated as dF = —2tt dn^ = 2nk. 



For a conserved electric current, 5j — 0, we call 



W A \j]=exp[i(6,j)] 



(2.3) 
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the Abelian Wilson loop. Its electric-photon (Wph) and the magnetic-monopole (Wmo) parts are specified as follows. 
Applying the Hodge decomposition to 9, we have the relation 

= A _1 A0 

= A- 1 (Sd + dd)9 

= A-^F + 2ttA- W m ) + A- 1 ^ 

= A^SF + 2TrA- 1 Sp + 27rq + A- 1 dS(9-27rq). (2.4) 

In the last line, we have used the relation rv- m ' — P + dq, where p{C2), q{C\) <E /Z. This means that an arbitrary 
shape of the open magnetic Dirac string rS m ^ is in general described by the sum of a fixed open string p with dp = —k 
and the closed strings dq with d 2 q — 0. Since all possible closed string fluctuations are summed over, one can choose 
an arbitrary open string p. Inserting Eq. I|2.4|l into Eq. 1)2. 3J1 , the Abelian Wilson loop can be written as 

W A \j) = exp^A- 1 ^, j)} ■ exp[ l (2^A- 1 ^,j)] = W Ph \j] ■ W Mo [j}, (2.5) 

where the third and fourth terms of Eq. (|2.4ll do not contribute to this decomposition because of the relations 
exp{2iri(q, j)] = 1 and 5j = 0. 

Let us proceed with the path integration of the partition function l|2.1|l keeping track of the two parts of the Wilson 
loop, Wp/Jj] and Wmo[j]- For simplicity and for picking up the essence of the following discussions, we restrict the 
differential operator in Eq. I|2.1|l to the leading term, D = /3 e A _1 . The path integral duality transformation of such 
a model itself has been discussed in many places since the works 

We first rewrite the summation over Dirac strings as the independent summation over monopole currents k (with 
constraint dk = 0) and q as 

E = E E- (2-6) 

n< m )e^ fee^, dk=o qez: 
Then, the integration with respect to 9 is replaced by 

m E = / mPh > ( 2 - 7 ) 

where 9 ph = A^ 1 SF represent noncompact link variables. Acting with an exterior derivative on 9 mo = 2irA~ 1 5p, one 
finds d9 mo = 2n(n^ + C (m) - dq) with C (m) = A^Sk. Thus, the partition function is written as noncompact QED 
with summation over closed monopole currents, 



Z = 



[ V9 ph exp -^.(d9 ph + 2ttC^) 2 +i{9 ph + 9 mo ,j) 



(2.8) 



In this expression one still realizes the violation of Abelian Bianchi identity in the form dF = 2ndC i - m ^ = 2nk 
due to dC<» = k. Using the relation (d9 ph : C^) = {9 ph ,5C^) = (since 8C {m ^ = 0) one can write (F) 2 = 
(d9 ph ) 2 + 47r 2 (C'"^) 2 . Taking into account the gauge fixing condition 59 ph = 0, one can integrate over 9 ph . This 
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yields a direct interaction term between electric currents j via the Coulomb propagator A 1 . Thus we have 



H exp 

ke%, dk=0 



(2.9) 



Defining C^(*C2) = A 1 S * j in analogy to = A 1 8k, the first term of the action can also be written in the 

form 



(2.10) 



where we have introduced the (inverse) dual gauge coupling (3 m = 1/g 2 , which should satisfy 47r 2 /3 e /3 m = 1 (i.e. 
Dirac's condition eg = Ait). Similarly, the square of C^ m ' can be rewritten 



^ft(C< m ») 2 = -±-(k,A^k) = -^(*k,A-Uk) 



(2.11) 



The exponential of this expression can further be understood as resulting from functional integration over the magnetic 
part of a noncompact dual gauge field 9 mo (*Ci), minimally coupled to the magnetic monopolc current, 



exp 



1 1 r°° 

— (♦fc.A- 1 **;) = / V9 mo cxp 

2Pm J — ao 



-y(dr) 2 + i(r,*i:) 



(2.12) 



We have attached the superscript " mo " in order to distinguish it from the photon part of the dual gauge field, 
9 ph (*Ci), which is defined in analogy to 9 mo as 



9 ph = 2 7 rA- 1 (5n (e) . 



(2.13) 



Here (*C2) G -ZT denotes electric Dirac strings, satisfying dn^ = which necessarily accompanies the presence 
of external electric charges. Formally, 9 ph enters our consideration when we re-express the monopole Wilson loop, 
using the relation 



(6> mo ,j) - (9 ph ,*k) = -2ir(p,*n^) = 2irN (N e Z) 



(2.14) 



This means that the direct coupling of j to 9 mo can be set equal to that of k to 9 ph , because of exp[i(27r7V)] = 1. 
Thus, the partition function is found to be 



f 

J — ( 



V9 mo ex P 
ke%, dk=a 



(3n 



id9 mo + 2irC {e) ) 2 + i{9 mo + 6 ph , *k) 



(2.15) 



The action is invariant under the transformation 9 mo \— » 9 mo + df. This is nothing but the realization of the dual gauge 
symmetry, due to the conservation of magnetic monopole currents, dk = 0. In this action, the electric currents are 
now implicitly defined via the violation of the dual Abelian Bianchi identity written down for the dual field strength 



F = d9 mo + 27rC (e) 



(2.16) 
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as dF = 2TtdC^ =2n*j, where dC^ 



The summation over monopole currents is the most difficult part of the evaluation. In principle, one needs to know 
the monopole dynamics, for instance, such as monopole current distribution in the vacuum and self interactions, etc. 
The numerical investigations of the effective monopole actions based on lattice Monte-Carlo simulation in the MAG 
provide such information, which has suggested the approximate form of the AP action given in Eq. (J2HJ. Here, we 
are not going to deal with these complication, since the kinetic structure of the dual gauge field, being composed 
of a regular 8 mo and a singular 8 ph parts, is not affected by the summation over monopoles. We then assume that 
the monopole current system is described by the grand canonical ensemble of closed loops, interacting via the dual 
gauge field. Then the complex-valued scalar monopole field Xi which minimally couples to the dual gauge field, is 
introduced 

HQ 

instead of monopole currents as 



£ exp 



i(9 mo + 8 ph ,*k) 



V\Vx exp 



{\(d + l (e mo + 9 ph )) x \ 2 + x(\x\ 2 -v 2 ) 2 } 



(2.17) 



where the A|x| 4 (A > 0) term plays the role to keep the density of loops being finite (it produces a short distance 
repulsion between the loop segments) and v denotes monopole condensate which describes the typical scale of the 
system. In this way, we arrive at the DAH model, 

Sdah = ^-{Ff + \(d + i0 m ° + ~8 ph )) x \ 2 + A(|x| 2 - v 2 ) 2 . (2.18) 

Although we cannot argue the precise values of the effective couplings in this formal derivation, we can restrict 
ourselves to the range of parameters able to describe the condensed phase of monopoles, according to the lattice 
results 0. 

Due to the singular structure of 8 ph associated with (see, Eq. (|2.13|) ). the DAH model has the open flux-tube 
solution, obtained by solving the field equations, 

PMdJ™° - dj™°) = 2(c° + 8 ph )^ 2 = k , (2.19) 

dM + (0™° + ~ d t)^ = 2A0(0 2 - v 2 ) . (2.20) 

Here, we have inserted the polar decomposition of the monopole field x = <fiexp(ir]) (</>, 77 £ 3?), and the phase r\ has 
been absorbed into the definition of The boundary conditions of the dual gauge field and monopole field are 

determined so as to make the energy of the system finite: just on the electric Dirac string rS e \ 8™° — and = 
whereas at large distance from the string, = —Off 1 and (j> — v. After solving the field equations (in general, 
numerically) , we can compute the profile of the electric field as the spatial part of the field strength F in Eq. (|2.16|l , 

E = V x e m0 + 2nC {e) = E mo + E ph , (2.21) 

and the magnetic current as the spatial part of the monopole current, 

~mo ~ph y o 

k = 2(0 +6 )cj) 2 , (2.22) 

respectively. Concrete forms of the field equations and the boundary conditions of fields for the straight q-q system 
are given in the Appendix 1X1 
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A typical flux-tube solution, the profile of the electric field and the monopole current, for the straight q-q system is 
shown in Fig.Q] The parameters we have chosen are (3 m — l/g 2 — 1, ra^-a — \/2gv-a — 0.5 and m x -a = 2y\va = 0.5, 
taking the q-q separation r = 16 a, where a is a certain length scale. The Ginzburg-Landau parameter is here 
k = m x /rriB = 1, which means the vacuum has superconducting properties just between type-I and type-II vacuum. 
This set of parameters is just to illustrate the flux-tube profile as an example. In Fig.[5J we then show the ingredient 
of the electric- field profile based on Eq. (|2.21(1 . We plot the flux-line pattern of the electric fields along the q-q axis, 
and the strength of each field as a function of the cylindrical radius. In Fig. we plot the strength of azimuthal 
monopole-current profile. Here, for the plots of the electric field and monopole current, we have added two cases 
corresponding to raj • a = 0.5 and m x ■ a = 0.25 (k — 0.5 : type I) and mg • a = 0.5 and m x ■ a = 0.75 (k = 1.5 : type 
II). Note that only the monopole-related part depends on k, while the photon part does not. 

In Fig. [21 we find that although the electric field derived from the monopole part of the dual gauge field, E™°, 
takes positive value near the center, it turns negative beyond a certain radius p c (in the given case, p c ~ 7a): this 
signals the appearance of a solenoidal electric field which plays an important role to cancel the Coulombic field, E ph , 
induced by electric charges, at some distance from the electric Dirac string. By this interplay the total electric field 
is squeezed from the dual superconducting vacuum, which finally leads to a flux tube. This is the composed internal 
structure of the DAH flux tube we are referring to. The solenoidal electric field and monopole supercurrent are related 
by the relation, V x E mo = k. It is important to realize that although the shape of total electric field profile becomes 
steeper as increasing k due to the change of its monopole part, the flux tube is always composed of the Coulombic 
and solenoidal electric fields. For the infinitely separated q-q system, the Coulombic contribution disappears and only 
the solenoidal electric field remains, where translational invariance of the flux-tube profile along the q-q axis becomes 
manifest. 

Now we come to the main point of the present section. Through the path-integral duality transformation, which 
has formally led us to the DAH model, we have found the role of the photon Wilson loop Wph[j] = exp[i(6 ph , j)] 
and the monopole Wilson loop Wmo[?] = exp[i(6 mo , j)] for the DAH model and its flux-tube solution; Wp/Jj] leads 
to the square of the Coulombic field strength after the integration over 8 ph (see, Eq. l|2.1U|l 1. while WmoL?'] i s 
translated into the interaction term between 6 ph — 2ir /S.~ x 5n^ e ' and the monopole field \ ( see i Eqs. (|2.14|) and 12.17|) '). 
Namely, the photon Wilson loop provides the origin of the Coulombic electric field contribution to the DAH flux tube. 
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FIG. 2: The flux-line pattern of the electric field (upper row) and the electric field strength as a function of the cylindrical 
radius (lower row): (a) the solenoidal electric field E mo and (b) the Coulombic field E ph add up to the flux-tube profile of the 
full electric field (c). 
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FIG. 3: The monopole current strength as a function of the cylindrical radius. 



On the other hand, the monopole Wilson loop determines the non-trivial behavior of the dual gauge field, inducing 
the monopole supercurrent and the solenoidal electric field component of the DAH flux tube. The Coulombic and 
solenoidal electric fields are responsible for the Coulombic and linearly rising parts of the inter-quark potential in the 
DAH model. In the actual AP lattice gauge simulations, it has been numerically shown that the potential detected 
by the photon and monopole Wilson loops have just the same feature Now, this is naturally understood 

from the relation between each Wilson loop and the composed internal structure of the DAH flux tube. We then 
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expect that the AP flux tube will exhibit the same composed structure as in the DAH flux tube, where Wp/jfj] and 
V^Mob] would be respective sources. 



III. DETECTING THE COMPOSED STRUCTURE OF ABELIAN-PROJECTED FLUX TUBE 



In this section, we are going to confirm the composed structure of the AP flux tube by measuring the electric field 
and monopole current profiles induced from the photon and monopole parts of the Abelian Wilson loop, based on the 
Monte-Carlo simulation of SU(2) lattice gauge theory in the MAG. 

In order to measure the the flux-tube profile induced by the Abelian Wilson loop W^aL?] = exp[i(#, j)], one can 
schematically use the following relation for a local operator O: 



Q^SnW O exp[-±(F,ADF)+t(6,j)} 
Z*™Y. n w exp[-i(F,ADF) + i{e,j)] 

S-„ m [J] cx P [- \ (F, ADF)] 

/_„ Z>0 E»c») W A [j] exp [- 1 (F, ADF)} 
= (OW a [j})q 
(W A \j])o ' 



where (• • -)j denotes an average in the vacuum with an external source, and (• • -)o an average in the vacuum without 
such source. Thus by measurement of the expectation values of (OWa)o and the Abelian Wilson loop (Wa)o, the 
expectation value of a local operator associated with the external source, (0)j, can be evaluated. Below, the Abelian 
field strength F and the monopole current k have been chosen as local operators O. In the actual simulation, since 
we do not know the exact form of the AP action, we first generate non-Abelian SU(2) gauge configurations and then 
specify the U(l) degrees of freedom by Abelian projection after MAG fixing. 

Typical profiles of the electric field and monopole current measured in this context are shown in Fig. 01 (some details 
of the simulation are given below briefly and in Appendix 0). Already at glance, the shape of the resulting profiles 
are very similar to the flux-tube profiles obtained within the DAH model, see Fig. ^ 

Before discussing the numerical simulation further, it is useful to consider what happens if we insert Wph and Wmo 
into Eq. [EHJl instead of W A - Using Eq. (|2~5|) and writing the Abelian field strength as F = d9 ph +2irC {m) = F Ph +F Mo , 
we can expect 



(JA ((F Ph + F M o)Wp h [j}WMo[j})o „ (Fp h W Ph [j}}o , (F Mo W M o[j})o .„ . , ,„ . _v 

~ 7777 — TTTTT r~T\ 7T77 — r~T\ 1 7777 FT\ = \*Ph)j + \^Mo)j ■ [o.Z) 

(W Ph [j]W M o[j})o {W Ph []])a {W Mo []])a 

Here, we have taken into account that in many cases lattice simulations in the MAG have found operators Xph and 
Ymoi defined in terms of the photon part and the monopole part of the Abelian link variable 9, respectively, to be 
uncorrelated: (XphYMo)o ~ (Xph)o(YMo)o (see, e.g., Ref. 40] and references therein). From the relation l|3.2f> . we 
expect that the sum of the flux profiles induced by the photon and monopole Wilson loops reproduces the total AP 
flux tube. 



Second, let us consider the expectation value of the monopole current k. Since we have the obvious relation 
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FIG. 5: Electric field profile from correlators with Abelian, photon and monopole Wilson loops at r — 6a — 0.48 fm (left) and 
at r = 12a = 0.97 fm (right). 



k = —dC( ,n ) = tcMo) where a photon part of the monopole current does not exist, kph oc d 2 9 ph = 0, we will observe 

/n _ (kMoW M o[j})a , , . 

( fc >J ~ — TTTr r~T\ = (kMo)j- (3.3) 

{vVMo[]\)o 

This means that the correlator of the monopole current only with the monopole Wilson loop will account for the 
full expectation value of monopole current profile and, at the same time, the correlator with the photon Wilson loop 
vanishes everywhere. 

We then show the corresponding lattice results, the electric field profile in Fig. [S] and the monopole current profile 
in Fig. El both as a function of the cylindrical radius. These measurements have been done at (3 = 2.5115 on a 32 4 
lattice after the MAG has been fixed. The q-q distances are r = 6a and 12a, and the measurements refer to the 
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FIG. 7: The same plots as in Fig. for flux tubes of length r = 0.48 fm and r = 0.97 fm, with the E z axis rescaled. The 
profile directly from the Abelian Wilson loop is omitted. 



x-y plane at half-distance. The lattice spacing is a = 0.081 fan, which has been determined from the non-Abclian 
string tension er^, ^Ja p hy S = y/^T/a = 440 MeV. Physically, the q-q distances correspond to 0.48 fm and 0.97 fm, 
respectively (see, Appendix iBl) . In Fig. we show the same electric held profiles as in Fig. EI focussing on the region 
where the monopole part of E z becomes negative. 

We hnd that these lattice results concerning the behavior of the profiles strongly support our considerations above; 
from the photon and the monopole Wilson loops, we obtain the Coulombic electric held and the solenoidal electric 
held with the monopole supercurrent profile, respectively. We find that the sum of these two contribution reproduces 
the profile obtained from the complete Abelian Wilson loop (see, Eq. (|3.2(l ). There is no correlation between the 
photon Wilson loop and monopole current as anticipated in Eq. (|3.3J) . Hence, we conclude that the AP flux tube has 
the same composed structure as the DAH flux tube. 

The behavior of the profiles as a function of the q-q distance r is also remarkable. While the monopole Wilson loop 
contributions, the solenoidal electric field and the monopole current profiles in the midplane arc rather stable with 
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respect to r, the photon Wilson loop contribution (i.e. the Coulombic electric field) drastically changes. From Fig. [5] 
it becomes obvious that the latter determines the change of the full Abelian electric field profile for different r. In 
order to see a really translationally invariant profile of the electric field, we need practically infinite q-q separation, 
r — > oo. In this limit, the profile only from the monopole part remains. This situation is also the same in the DAH 
flux tube. 



IV. SUMMARY 

It has already been known that the profiles of the classical flux-tube solution in the dual Abelian Higgs (DAH) 
model and of the Abelian-projected (AP) flux tube, observed in lattice simulations in the maximally Abelian gauge 
(MAG), look quite similar. 

In this paper, in order to establish a more detailed correspondence between these two kinds of profiles, we have 
studied the composed structure of both flux tubes more carefully. First, by applying the path-integral duality trans- 
formation to the Villain type compact QED considered as the approximate action of the AP gauge theory, we have 
been led to the U(l) DAH model. Along the way, we have identified the electric and magnetic parts of the Abelian 
Wilson loop by the Hodge decomposition, and have clarified the role of each contribution to the structure of the 
flux-tube solution in the DAH model. The photon and monopole Wilson loops provide sources of the Coulombic and 
solenoidal electric field components of the DAH flux tube. 

Guided by this observation, we have performed lattice simulations of the SU(2) lattice gauge theory in the MAG 
and have measured the flux profiles induced by the photon and the monopole Wilson loops. We have found that the 
resulting profiles show the same composed structure as the DAH flux tube. 

The further question would be how both sides are related quantitatively. One way would be to fit the profile of 
the AP flux tube by that of the DAH flux tube and to determine the DAH parameters which remain unknown in 
the formal derivation of the DAH model. Here, we would like to emphasize that the composed structure of the AP 
flux tube found here and its relation to the DAH flux tube will be important for further quantitative discussions. In 
fact, there is no such a work that takes into account the correspondence of the structures. In addition to this, as 
we have mentioned briefly in the introduction, a more systematic study of the AP flux-tube profile itself is required: 
the Gribov copy effect in the MAG, the scaling property, the q-q distance dependence, etc. Otherwise, one cannnot 
trust the robustness and p hysi cal relevance of the resulting DAH parameters. A par t of such a quantitative analysis 
is reported in Lattice 2002 |46j and the detailed report will be presented in Ref. 41] . 

In closing, we note that although we have concentrated here on SU(2) gauge theory, the ideas discussed in the 
present paper can be extended to arbitrary AP-SU(iV) gauge theory in the MAG 
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APPENDIX A: CLASSICAL FLUX-TUBE SOLUTION IN THE DAH MODEL 



In this appendix, we present concrete form of the field equations of the DAH model, Eqs. I|2.19[) and (|2.20() . for a 
straight q-q system. We use here the continuum notations; let B™° be the continuum form of the regular dual gauge 
field denoted 0™° on the lattice, Bf that of the singular dual gauge field Of/ 1 . 

We put the quark and the antiquark at x\ = (—r/2)e z and = (+r/2)e z . The corresponding electric current 
is written as j^x) = <5 M o (S(x — x\) — S(x — a^))- Since this system has cylindrical geometry, the fields can be 
parametrized in terms of cylindrical coordinates (p, tp, z) as 



4> = <j){p,z) 



B ph 



B m °{p,z)e, 



_ B mo (p,z) 



z + r/2 



-tp ) 



z - r/2 



2p V y/(P + (z + r/2) 2 VP 2 +¥ Z r/2) 2 



(Al) 
(A2) 

(A3) 



The factor n in B ph is the win ding number of the flux tube (an integer value) , which is determined by the representation 
of the electric charges |48|, |50j, |51| . The fundamental representation corresponds to n = 1. 

The field equations l|2.19|l and 12.20|) are then reduced to: 



0T, 



d 2 B r ' 



1 dB mo d 2 B mo 



2 B 



dp 2 p dp dz 2 
z + r/2 



z - r/2 



^p 2 + (z + r/2) 2 ^p 2 + {z-r/2f 

d 2 4> 1 d<t> d 2 (j) 
dp 2 p dp dz 2 

( jjmo _ n ( z+r/2 _ z-r/2 \ \ 




j^mo n / z-t-r / z z — r/2 



= o, 



<t>- 2A0(0 2 - v z ) = 0. 



(A4) 



(A5) 
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The boundary conditions are specified so as to make the energy of the system finite as 

B mo = as p -> 0, and = as p -> for - r < z < r, 
z + r/2 z-r/2 
~ \ \jp 2 + (z + r/2) 2 yfp? + (z-r/2) 2 



B™ = 2 ( Z + r/2 - Z r/2 I and ,,= ,• as ; -m . (A(,l 



After getting the numerical solution of the field equations for B mo and 0, the profiles of the electric field are computed 
as Eq. I|2.21[) . where 

ldB mo ldB mo 

E = ~~ — aZ — e P + ~ — ST - ez > ( A7 ) 



E P 



p dz p dp 

P h _ n ( P_ 



2 \(p 2 + (z + r/2) 2 ) 3 / 2 (p 2 + (z- r/2) 2 ) 3 / 2 
n ( z + r/2 z — r/2 



2 \(p 2 + (z + r/2) 2 ) 3 / 2 (p 2 + (z- r/2) 2 ) 3 / 2 
The profile of the monopole current l|2.22J) is given by 



(A8) 



* = 2 | U - J [ Z + r/2 2 - — /2 2 || . 

Yp3 + (« + r/2) 2 + (z - r/2) 2 



APPENDIX B: LATTICE SIMULATION DETAIL 



For the SU(2) link variables, U^s), generated by Monte-Carlo method with Wilson gauge action, we adopt the 
maximally Abelian gauge (MAG) fixing, which is achieved by maximizing the functional 

R[U v }=J2^{r 3 Uj(s)r 3 U^ t( s )} . (Bl) 

After the MAG fixing, Abelian projection is performed; the SU(2) link variables U^(s) = U^f A (s) are factorized 
into a diagonal (Abelian) link variable w M (s) <E U(l) 3 and the off-diagonal (charged matter field) parts c M (s), e*(s) <E 
SU(2)/U(1) 3 as follows 



(B2) 



TT MAf, _ ( V 1 - M*)l 2 - c U s ) \ ( o \ 

where the Abelian link variables u^s) are then explicitly written as 

«„(*)= e""« (^(«)e[-7r,7r)). (B3) 
The Abelian plaquette variables is then constructed as 

0„„(«) = e M (a)+6 v {a + lj.)-e ii {a + 0)-0 v (a) e [-An, 4*) , (B4) 
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which is decomposed into a regular part 0^„(s) S [— n, n) and a singular (magnetic Dirac string) part n^\s) = 
0,±1,±2 as follows 

M s ) = 5 lw (s)+2wn$>(s) . (B5) 

The Abelian field strength is defined by 9^(s) = Q^ v (s) — 2Ttnffi (s). Following DeGrand and Touissaint |S2|, magnetic 
monopoles are extracted as the string boundaries 

M s d) = - \ e ^P°dv n{ °po ) ( s + A) (ei234 = 1) , (B6) 

where |fc M (sd)| < 2 and Sd = s + (l + 2 + 3 + 4)/2 denotes the dual site. 

For measuring the correlation function, we have used the following local operators: an electric field operator 

Q(s) = i9 i4 (s) = i(9 i4 (s) - 2^4™ ) ( s )), (B7) 
and a monopole current operator 

0{s d ) = 2niki(sd) (B8) 
The Abelian Wilson loop is constructed as 

W^fol = ]!"/*(*) = e iE,<yM * ) - (B9) 
lej 

Similarly, the photon and the monopole Wilson loop are constructed from the photon and monopole parts of Abelian 
link variables, 8 ph and 9 mo , respectively, where 

6,( S ) = A- l d u (9„ v (s) + 27m#>00) = Ofis) + 9™(s). (BIO) 

In this decomposition, it is necessary to adopt the Abelian Landau gauge which is characterized by d^9^(s) = 0. 
Note, however, that the Wilson loops constructed from each link variable are Abelian gauge invariant. 

In this simulation, in order to see the profiles which belong to the ground state of a flux tube, we have adopted 
a smearing technique for spacelike Abelian link variables. Then we have constructed the smeared Abelian Wilson 
loop • Considering the fourth direction as the Euclidean time direction, we have performed N s times the following 
step in a smearing procedure applied only to the spatial Abelian links (i,j — 1, 2, 3), 

ae i6i{s) + y^ e i(0i(a)+0i(a+i)-0i(«+i)) gifcOO ) ( BU ) 

where a is an appropriate smearing parameter. The same procedure was also applied to the spatial parts of the 
photon and the monopole link variables before constructing each type of Wilson loop. 
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The numerical simulations which are presented in this paper have been done at {3 = 2.5115. The lattice volume was 
32 4 . We have used 100 configurations for measurements. We have produced them after 3000 thermalization sweeps, 
separated by 500 Monte Carlo updates. They have been stored for performing MAG fixing. This has been repeated 
N g times, starting each time from a different random gauge copy of the configuration, in order to explore an increasing 
number of Gribov copies. The copy reaching the maximal value of the gauge functional (|Blfl has been selected for 
measuring the profiles and kept for further increasing of N g . Finally we have chosen N g = 20. For the MAG fixing 
itself, we have used the simulated annealing algorithm 0], followed by a final steepest descent relaxation. The size 
of the Wilson loops mainly studied (for Fig. EJ El and 0) are R x T = 6 x 6 and 12 x 6 in units of lattice spacing a. 
We have measured the profiles in the x-y plane orthogonal to the Wilson loop in its midpoint. The Abelian smearing 
parameters have been found by optimization as N s = 8 and a = 2.0. With this choice, the profiles induced by the 
Abelian Wilson loop with timelike extensions T = 8 and T — 6 agree within errors. 

The physical scale (the lattice spacing a(/3 — 2.5115)) has been determined from the non-Abelian string tension 
o~ l by fixing ^/a p h ys = y/o'L/a = 440 MeV. The non-Abelian string tension has been re-evaluated by measuring 
expectation values of non-Abelian Wilson loops with an optimized non-Abelian smearing. The potential has been 
fitted to match the form V(R) = C — A/R + ctl R. The resulting string tension is (Tl = 0.0323(4) at (3 = 2.5115, such 
that the corresponding lattice spacing in physical units is a(f3) — 0.0806(5) fm. 
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